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PREFACE 


The  screening  analysis  reported  herein  is  based  on  seismic  stability 
evaluations  of  several  earth  dams,  in  particular  Richard  B.  Russell  Dam  in 
Georgia  and  South  Carolina  and  Ririe  Dam  in  Idaho,  and  was  performed  by  the 
Earthquake  Engineering  and  Geophysics  Division  (EE&GD),  Geotechnical  Labora¬ 
tory  (GL),  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES),  Vicksburg, 
Miss.  This  study  was  sponsored  by  the  Office,  Chief  of  Engineers  (OCE),  U.  S. 
Army,  under  the  Civil  Works  Investigational  Studies  (CWIS),  Soils  Research 
Program  (Work  Unit  31145),  "Liquefaction  of  Dams  and  Foundations  During  Earth¬ 
quakes,”  for  which  Mr.  R.  F.  Davidson  was  the  OCE  Technical  Monitor. 

The  research  was  conducted  and  the  report  prepared  by  Ms.  M.  E. 
Hynes-Griffin,  EE&GD,  and  Dr.  A.  G.  Franklin,  Principal  Investigator  and 
Chief,  EE&GD.  Appendix  A  was  prepared  by  Mr.  F.  K.  Chang,  EE&GD.  The  study 
was  performed  under  the  general  supervision  of  Dr.  W.  F.  Marcuson  III,  Chief, 
GL. 

COL  Tilford  C.  Creel,  CE,  was  Commander  and  Director  of  WES  during  the 
preparation  of  this  report.  Mr.  F.  R.  Brown  was  Technical  Director. 
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RATIONALIZING  THE  SEISMIC  COEFFICIENT  METHOD 


PART  I :  INTRODUCTION 

1.  Until  the  1960's,  seismic  analysis  of  dams  consisted  essentially  of 
the  seismic  coefficient  method,  in  which  a  static,  horizontal  inertia  force 
was  applied  to  the  potential  sliding  mass  in  an  otherwise  conventional  static 
limit  analysis.  The  magnitude  of  the  inertia  force  was  chosen  on  the  basis  of 
judgment  and  tradition;  a  rational  basis  was  lacking.  Alternatives  to  this 
approach  became  available  during  the  1960's  and  1970' s.  A  method  of  analysis 
that  dealt  with  the  softening  or  liquefaction  of  granular  soil.?  was  evolved, 
largely  through  work  at  the  University  of  California  at  Berkeley,  with  Profes¬ 
sor  H.  B.  Seed  playing  the  leading  role.  This  approach  is  based  on  comparison 
of  dynamic  shear  stresses  computed  in  a  transient  response  analysis  to  the 
cyclic  strength  (resistance  to  liquefaction)  obtained  from  laboratory  cyclic 
shear  tests  or  from  empirical  correlations  of  liquefaction  occurrence  (and  non¬ 
occurrence)  with  Standard  Penetration  Tests  (Seed,  et  al.  1975a,  b;  Seed  1979; 
Seed  and  Idriss  1983).  A  second  alternative  is  to  deal  with  the  permanent  de¬ 
formations  that  might  be  anticipated  if  the  embankment  and  foundation  soils  do 
not  suffer  liquefaction  or  severe  softening  under  cyclic  loading,  using  as  an 
idealized  model  of  the  displaced  part  of  the  embankment  a  rigid  block  sliding 
on  an  inclined  plane.  This  approach  was  proposed  by  the  late  Professor  N.  M. 
Newmark  in  his  Rankine  Lecture  (1965).  Other  contributions  to  a  coherent  pro¬ 
cedure  using  this  approach  have  been  made  by  Professors  Ambraseys  and  Sarma, 
Imperial  College,  London  (e.g.  Ambraseys  and  Sarma  1967;  Sarma  1975,  1979), 

the  Berkeley  group  (Goodman  and  Seed  1966,  Makdisi  and  Seed  1977),  and  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  (Franklin  and  Chang 
1977,  Franklin  and  Hynes-Griffin  1981). 

2.  Sufficient  experience  has  been  gained  in  the  application  of  the  New¬ 
mark  approach  to  allow  some  conclusions  to  be  drawn.  In  the  absence  of  lique¬ 
faction  effects,  dams  with  adequate  static  factors  of  safety  against  sliding 
are  not  likely  to  be  predicted  by  this  analysis  to  be  subject  to  deformations 
so  large  as  to  endanger  their  reservoirs,  through  limited  sliding  deformations 
may  be  predicted.  This  result  suggests  that  many — perhaps  most--perroanent 
displacement  analyses  do  not  really  need  to  be  done,  and  that  some  simple 
screening  method  should  be  applied  to  separate  those  dams  that  are  clearly 


safe  against  earthquake-induced  failure  from  those  that  require  further  anal¬ 
ysis.  A  seismic  coefficient  analysis  can  serve  this  screening  function,  be¬ 
cause  the  accumulated  experience  in  permanent  displacement  analyses  now  pro¬ 
vides  a  rational  basis  for  choosing  the  value  of  the  coefficient.  The 
rationale  is  based  on  assuring  that  deformations  will  be  limited  to  tolerable 
values,  assuming  the  worst  combination  of  earthquake  loads  and  resonant  em¬ 
bankment  response.  A  procedure  that  uses  this  approach  is  proposed  in  this 


PART  II:  PERMANENT  DISPLACEMENT  ANALYSIS 

3.  The  major  components  of  a  permanent  displacement  analysis  of  the 
Newmark  type,  as  applied  by  the  WES,  are  shown  in  Figure  1.  The  primary  com¬ 
ponent  is  the  analysis  of  motions  of  a  system  consisting  of  a  rigid  block 
sliding  on  an  inclined  plane,  chosen  to  represent  a  potential  sliding  mass  in 
an  embankment,  as  described  by  Newmark.  A  conventional  limit  analysis,  or 
slope  stability  analysis,  with  slight  modifications,  provides  the  shearing 
resistance  between  the  block  and  plane.  Because  bedrock  motions  may  be  ampli¬ 
fied  upon  being  propagated  upward  through  an  embankment,  a  rigid-body  model 
may  underestimate  displacements,  and  an  analysis  of  the  amplification  response 
of  the  embankment  is  incorporated  to  account  for  amplified  accelerations  in 
the  embankment. 


Figure  1.  Permanent  displacement  analysis 


Stability  Analysis 

4.  The  concept  of  the  traditional  pseudostatic,  seismic  coefficient 
method  of  analysis  is  illustrated  in  Figure  2.  In  an  otherwise  conventional 
static  stability  analysis,  such  as  a  method  of  slices  analysis,  the  earthquake 
loading  is  represented  by  a  statically  applied  horizontal  force  kW  ,  where  W 
is  the  weight  of  the  slice  and  k  is  the  seismic  coefficient,  which  is  some 
fraction  of  gravity.  The  value  of  k  is  generally  prescribed  by  code  or  reg¬ 
ulation,  with  values  usually  in  the  range  of  0.05  to  0.20,  depending  on  the 
seismicity  of  the  site.  The  procedure  is  described  in  EM  1110-2-1902  (U.  S. 
Army,  Office,  Chief  of  Engineers  1970)  and  in  many  standard  texts. 


.  • 


.  • 
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Figure  2.  Earthquake  stability  analysis  by  pseudostatic 
method  using  seismic  coefficient 

5.  For  analysis  of  permanent  displacements,  the  shearing  resistance  be¬ 
tween  the  potential  sliding  mass  and  the  underlying  base  is  evaluated  in  terms 
of  a  critical  acceleration  N  ,  defined  as  the  acceleration  (of  the  ground  or 
embankment  below  the  sliding  surface)  that  will  reduce  the  factor  of  safety 
against  sliding  to  unity,  i.e.,  that  will  make  sliding  imminent.  The  value  of 
N  ,  which  is  expressed  as  a  fraction  of  gravity  (g)  ,  is  obtained  through  a 
stability  analysis  similar  to  conventional  pseudostatic  stability  analyses, 
but  which  includes  two  special  features.  One  is  that  the  stability  is  evalu¬ 
ated  in  terms  of  a  critical  acceleration  rather  that  a  factor  of  safety,  and 
the  other  is  that,  because  the  amplified  accelerations  vary  over  the  height  of 
the  embankment,  critical  accelerations  are  determined  for  possible  sliding 
masses  whose  bases  lie  at  various  elevations  in  the  section  (Figure  3). 

6.  The  analysis  may  be  performed  using  conventional  stability  analysis 


.  • 
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Figure  3.  Critical  acceleration  as  a  function  of  elevation 

methods  such  as  those  of  Bishop  (1955)  or  Morgenstern  and  Price  (1965).  Trial 
values  of  acceleration  may  be  used  to  find  the  value  that  reduces  the  factor 
of  safety  to  unity.  The  Sarma  method  (Sanaa  1975),  which  employs  a  slip  sur¬ 
face  of  arbitrary  shape,  determines  the  value  of  N  directly. 

7.  In  principle,  the  analysis  can  be  performed  on  either  a  total  or  an 
effective  stress  basis,  but  the  problems  of  estimating  pore  pressures  induced 
by  cyclic  shearing  are  avoided  by  using  a  total  stress  analysis.  The  usual 
Corps  of  Engineers  practice  for  static  stability  analyses  is  to  use  a  compos¬ 
ite  shear  strength  envelope  based  on  the  S  test  (consolidated-drained)  at  low 
confining  pressures  and  the  R  test  (consolidated-undrained)  at  high  confining 
pressures  (Figure  4).  This  strength  envelope,  which  conservatively  takes  into 
account  possible  dissipation  of  shear-induced  negative  pore  pressures  that 
might  occur  in  the  field  but  cannot  occur  in  an  undrained  test  in  the  labora¬ 
tory,  is  recommended  for  pervious  soils.  For  soils  of  low  permeability,  in 
which  undrained  conditions  are  more  likely  to  exist  during  an  earthquake,  an 
undrained  (R)  strength  envelope  would  be  appropriate. 

8.  Makdisi  and  Seed  (1977)  point  out  that  substantial  permanent  strains 
may  be  produced  by  cyclic  loading  of  soils  to  stresses  near  the  yield  stress, 
while  essentially  elastic  behavior  is  observed  under  many  (>100)  cycles  of 
loading  at  80  percent  of  the  undrained  strength.  They  recommend  the  use  of 

80  percent  of  the  undrained  strength  as  the  "dynamic  yield  strength"  for  soils 
that  exhibit  small  increases  in  pore  pressure  during  cyclic  loading,  such  as 
clayey  materials,  dry  or  partially  saturated  cohesionless  soils,  or  very  dense 
saturated  cohesionless  materials. 


normal  stress  a 

Figure  4.  Composite  "S-R"  strength  envelope 


Sliding  Block  Analysis 


9.  Figure  5  presents  the  elements  of  the  sliding  block  analysis  (Frank¬ 
lin  and  Chang  1977).  The  potential  sliding  mass  in  Figure  5a  is  in  a  condi¬ 
tion  of  impending  failure,  so  that  the  factor  of  safety  equals  unity.  This 
condition  is  caused  by  the  acceleration  of  both  the  base  and  the  mass  toward 
the  left  of  the  sketch  with  an  acceleration  of  Ng  .  The  acceleration  of  the 
mass  is  limited  to  this  value  by  the  limit  of  the  shear  stresses  that  can  be 
exerted  across  the  contact,  so  that  if  the  base  acceleration  were  to  increase, 
the  mass  would  move  downhill  relative  to  the  base.  By  D'Alembert's  principle, 
the  limiting  acceleration  is  represented  by  an  inertia  force  NW  applied 
pseudostatically  to  the  mass  in  a  direction  opposite  to  the  acceleration. 

10.  Figure  5b  shows  the  force  polygon  for  this  situation.  The  angle  of 
inclination  0  of  the  inertia  force  may  be  found  as  the  angle  that  is  most 
critical;  this  is,  the  angle  that  minimizes  N  .  Its  value  is  usually  within 
a  few  degrees  of  zero,  and  since  the  results  of  the  analysis  are  not  sensitive 
to  it,  it  can  generally  be  ignored.  The  angle  (3  is  the  direction  of  the  re¬ 
sultant  S  of  the  shear  stresses  on  the  interface  and  is  determined  in  the 
course  of  the  stability  analysis.  The  same  force  polygon  applies  to  the  model 
of  a  sliding  block  on  a  plane  inclined  at  an  angle  p  to  the  horizontal  (Fig¬ 
ure  5c).  Hence,  the  sliding  block  model  is  used  to  represent  the  sliding  mass 
in  an  embankment. 
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a.  POTENTIAL  SLIDING  MASS 


b.  FORCE  POLYGON 
FOR  F.S.=1.0 


c.  SLIDING  BLOCK  MODEL  d.  FORCE-DISPLACEMENT  RELATION 


e.  COMPUTATION  OF  DISPLACEMENT 

Figure  5.  Elements  of  the  sliding  block  analysis 


11.  The  force-displacement  relation  in  Figure  5d  is  assumed  to  apply  to 
this  system.  The  force  in  this  diagram  is  the  inertia  force  corresponding  to 
the  instantaneous  acceleration  of  the  block,  and  the  displacement  is  the  slid¬ 
ing  displacement  of  the  block  relative  to  the  base.  It  is  usually  assumed 
that  resistance  to  uphill  sliding  is  large  enough  that  all  displacements  are 
downhill.  This  assumption,  in  addition  to  simplifying  the  calculations,  is 
both  realistic  and  conservative  (Franklin  and  Chang  1977). 

12.  If  the  base  (i.e.,  the  inclined  plane)  is  subjected  to  some  se¬ 
quence  of  acceleration  pulses  (the  earthquake)  large  enough  to  induce  sliding 
of  the  block,  the  result  will  be  that  after  the  motion  has  abated,  the  block 
will  come  to  rest  at  some  displaced  position  down  the  slope.  The  amount  of 
permanent  displacement,  which  will  be  called  u  ,  can  be  computed  by  using 
Newton's  second  law  of  motion,  F  =  ma  ,  to  write  the  equation  of  motion  for 
the  sliding  block  relative  to  the  base,  and  then  numerically  or  graphically 
integrating  (twice)  to  obtain  the  resultant  displacement.  During  the  time 
intervals  when  relative  motion  is  occurring,  the  acceleration  of  the  block 
relative  to  the  base  is  given  by 

_  _  /  vr\  cos  (B  -  0  -  <b) 

U  arel  “  (abase  )  cos  <t> 

=  (‘base  -  ")  •  « 

where 

afe^  =  relative  acceleration  between  the  block  and  the  inclined  plane 

abase  =  acceleration  of  the  inclined  plane,  a  function  of  time 

N  =  critical  acceleration  level  at  which  sliding  begins 

0  =  direction  of  the  resultant  shear  force  and  displacement,  and  the 
inclination  of  the  plane 

0  =  direction  of  the  acceleration,  measured  from  the  horizontal 
<|>  =  friction  angle  between  the  block  and  the  plane 

The  acceleration  a^ase  is  the  earthquake  acceleration  acting  at  the  level  of 
the  sliding  mass  in  the  embankment.  It  is  assumed  to  be  equal  to  the  bedrock 
acceleration  multiplied  by  an  amplification  factor  that  accounts  for  the 
quasi-elastic  response  of  the  embankment. 

13.  The  permanent  displacement  is  determined  by  twice  integrating  the 

relative  accelerations  over  the  total  duration  of  the  earthquake  record.  It  is 
assumed  that  <)>  ,  0  ,  and  0  do  not  change  with  time;  thus,  the  coefficient 
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a  is  a  constant  and  is  not  involved  in  the  integration.  In  the  final  stage 
of  analysis,  the  result  of  the  integration  is  multiplied  by  the  coefficient 
a  ,  the  determination  of  which  requires  knowledge  of  the  embankment  properties 
and  the  results  of  the  pseudostatic  stability  analysis.  For  most  practical 
problems,  the  coefficient  a  differs  from  unity  by  less  than  15  percent  (Fig¬ 
ure  6).  For  the  purposes  of  this  report,  a  value  of  unity  will  be  assumed. 


4>,  degrees 

Figure  6.  Values  of  the  coefficient  a 

14.  The  integration  can  be  easily  visualized  on  a  plot  of  base  velocity 
versus  time,  obtained  by  a  single  integration  of  the  acceleration  record  (Fig¬ 
ure  5e).  Since  the  slope  of  the  velocity  curve  is  the  acceleration,  the 
limiting  acceleration  Ng  of  the  block  defines  .the  velocity  curve  for  the 
block  by  straight  lines  in  those  parts  of  the  plot  where  the  critical  accel¬ 
eration  has  been  exceeded  in  the  base.  The  area  between  the  curves  gives  the 
relative  displacement. 

15.  In  this  analysis,  the  characteristics  of  the  sliding  mass  in  the 
embankment  are  represented  only  by  the  critical  acceleration  N  and  the  am¬ 
plification  factor,  the  latter  being  simply  a  constant  multiplying  factor. 

The  permanent  displacement  u  for  a  particular  earthquake  record  can  be  de¬ 
termined  as  a  function  of  N/A  ,  where  Ag  is  the  peak  value  of  the  earth¬ 
quake  acceleration,  and  the  u  versus  N/A  curve  can  be  determined  from  the 
earthquake  record  without  reference  to  a  particular  embankment. 
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16.  Kutter  (1982)  has  done  limited  experimental  testing  of  this  method 
by  means  of  model  embankments  shaken  by  simulated  earthquakes  in  the  Cambridge 
University  geotechnical  centrifuge.  For  these  tests,  the  sliding  block  model 
gave  poor  predictions  of  very  small  displacements  (<1  cm),  but  if  strength 
degradation  was  provided  for,  it  produced  good  predictions  when  the  displace¬ 
ments  at  prototype  scale  were  greater  than  about  1  cm. 

17.  The  following  example,  drawn  from  Kutter 's  test  results  for  embank¬ 
ment  model  D  and  earthquake  I  ,  demonstrates  the  application  of  the  dis¬ 
placement  calculation  procedures  described  herein.  If  the  yield  acceleration 
for  the  embankment  model  is  calculated  on  the  basis  of  80  percent  of  the  mea¬ 
sured  shear  strength,  and  the  measured  amplification  of  the  base  motion  is  in¬ 
cluded,  the  predicted  displacement  is  15.0  cm  and  the  corresponding  measured 
displacement  is  16.4  cm  at  the  prototype  scale. 

18.  Sliding  block  analyses  have  been  done  at  the  WES  for  348  horizontal 
earthquake  components  and  6  synthetic  records.  These  calculations  are  tabu¬ 
lated  in  Appendix  A.  The  results  are  summarized  in  Figure  7,  which  shows  the 
mean,  mean  plus  one  standard  deviation  (o),  and  upper  bound  curves,  for  all 
natural  records  and  all  synthetic  records  representing  magnitudes  smaller  than 
8.0.  (Caution  is  recommended  in  interpreting  these  curves  quantitatively  in 
terms  of  relative  probability,  because  the  data  base  is  biased  by  over¬ 
representation  of  a  single  earthquake,  the  1971  San  Fernando  earthquake,  which 
produced  records  at  many  locations.) 

19.  The  question  of  how  much  deformation  is  tolerable  has  no  single 
answer;  it  depends  on  such  factors  as  the  size  and  geometry  of  the  dam,  the 
zonation,  the  location  of  the  sliding  surface,  and  the  amount  of  freeboard 
available.  The  authors  have  arbitratily  chosen  1  m  of  permanent  displacement 
as  a  tolerable  upper  limit.  Such  a  deformation  would  surely  be  considered 
serious  damage,  but  it  could  be  tolerated  in  most  dams  without  immediately 
threatening  the  integrity  of  the  reservoir.  The  unusual  cases  where  a  dam 
could  not  tolerate  1  m  of  displacement,  because  of  small  freeboard  or  vulner¬ 
ability  of  critical  design  features  to  small  displacements,  should  be  evalu¬ 
ated  by  other  methods. 

20.  When  Figure  7  is  entered  at  100  cm  (1  m)  of  displacement,  the  cor¬ 
responding  N/A  value  is  0.17.  Thus,  deformations  will  be  limited  to  less 
than  1  m  of  displacement  if  the  critical  acceleration  is  as  much  as  0.17  times 
the  peak  acceleration  on  the  sliding  surface. 
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N/A  • 

Figure  7.  Permanent  displacement  u  versus  N/A  ,  based  on  348 
horizontal  components  and  6  synthetic  accelerograms 


Embankment  Response  Analysis 


21.  Amplification  of  ground  motions  in  the  embankments  may  be  examined 
by  analysis  of  a  shear-beam  model  of  the  embankment- foundation  system.  A 
closed-form  solution  has  been  obtained  by  Sarma  (1979)  for  the  problem  illus¬ 
trated  by  Figure  8.  The  model  considered  is  an  untruncated  triangular  wedge 
of  height  h^  with  a  shear-wave  velocity  and  density  ,  underlain  by 

a  foundation  layer  with  thickness  ,  shear-wave  velocity  ,  and  density 
p2  .  Both  the  wedge  and  foundation  are  linearly  visco-elastic  and  have  the 
same  damping  ratio  D  .  The  earthquake  motions  are  considered  to  be  rigid-body 
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Figure  8.  Mathematical  model  for  viscoelastic  shear 
beam  analysis  of  embankment  and  foundation  response 
by  the  Sarma  method 


motions  in  the  rock  underlying  the  foundation  layer,  and  it  is  assumed  that 
all  motions  are  horizontal  (hence,  a  shear-beam  model).  Shear-wave  veloc¬ 
ities  and  damping  values  are  chosen  so  as  to  be  consistent  with  expected 
strain  levels.  The  computation  of  accelerations  is  carried  out  in  the  time 
domain. 

22.  Geometry  and  material  properties  are  described  in  terms  of  the 
dimensionless  parameters  m  and  q  ,  which  are  defined  as 


P1S1 

m  =  — -g—  and 
P2b2 


Slh2 
'  "  S2hl 


(2) 


23.  For  use  with  the  sliding  block  analysis,  accelerations  are  averaged 
over  a  wedge  that  is  selected  to  be  approximately  equivalent  in  volume  and 
location  to  a  potential  sliding  mass  with  its  base  at  some  chosen  elevation, 
as  shown  in  Figure  9.  The  average  acceleration  acting  on  the  wedge  at  any 
instant  is  taken  as 


a 


av 


(3) 


where  a(z)  is  the  acceleration  of  the  area  element  dA  ,  at  elevation  z  , 
and  A  is  the  total  area  of  the  wedge. 

24.  The  largest  average  acceleration  that  acts  on  the  wedge  at  any  time 
during  the  earthquake  shaking  is  produced  as  the  output  of  the  computer  pro¬ 
gram,  and  the  ratio  of  that  acceleration  value  to  the  peak  bedrock  accelera¬ 
tion  is  taken  as  the  amplification  factor  for  the  wedge.  In  Figure  10,  values 
of  the  amplification  factor  are  plotted  against  the  embankment  fundamental  pe¬ 
riod  Tq  for  one  record  of  the  Parkfield  earthquake  of  27  June  1966.  Curves 
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Figure  9. 


o 

o  2.0 


ca 

o 

~  1.0 

a. 

E 

<o 


fundamental  period,  seconds 

Figure  10.  Amplification  curves  for  the  S  25  W  component, 

Temblor  No.  2  Record,  Parkfield  earthquake  of  27  June  1966 
(damping  =  20  percent) 

are  shown  for  wedges  with  their  bases  at  various  distances  y/hj  (defined 
in  Figure  9)  from  the  crest,  for  a  single  combination  of  m  and  q  values 
(m  =  0.5,  q  =  0.5). 

25.  Amplification  curves  have  been  obtained  from  27  strong-motion 
earthquake  records  and  a  wide  range  of  m  and  q  values  (representing  em¬ 
bankments  on  rock  and  on  foundation  layers  of  varied  thickness,  and  with  a 
variety  of  relative  embankment- foundation  stiffnesses).  Damping  values  used 
ranged  from  15  to  20  percent.  Also,  numerous  computed  amplification  values 
have  been  obtained  from  finite  element  analyses  and  from  the  literature. 
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Figure  11  presents  a  summary  of  computed  resonant  response,  obtained  by  plot¬ 
ting  the  values  at  the  peaks  of  the  amplification  curves.  Table  1  shows  these 
peak  amplification  values.  Amplification  values  obtained  from  finite  element 
analyses,  which  do  not  necessarily  represent  resonant  conditions,  are  gener¬ 
ally  lower  than  these  curves  indicate. 

26.  To  use  these  curves  in  a  permanent  displacement  analysis,  pick  off 
the  amplification  factor  for  the  depth  of  sliding  being  investigated,  and  mul¬ 
tiply  the  peak  bedrock  acceleration  by  that  value  before  entering  the  plot  of 
displacement  versus  N/A  .  This  step  involves  an  assumption  that  the  sliding 
block  analysis  and  the  amplification  analysis  can  be  decoupled.  In  fact, 
there  is  good  reason  to  believe  that  decoupling  results  in  overestimates  of 
the  amplification  when  very  strong  shaking  is  involved.  The  amplification  may 
be  large  in  cases  where  the  motions  are  small  and  the  embankment  behavior  is 
nearly  elastic  (Gazetas,  et  al.  1981),  but  this  assumption  is  not  compatible 
with  inelastic  embankment  response.  If  accelerations  are  high  enough  to  pro¬ 
duce  sliding  on  a  deep  surface,  then  the  embankment  is  incapable  of  propa¬ 
gating  these  large  accelerations  to  higher  elevations.  The  critical  accelera¬ 
tion  on  a  slip  surface  defines  the  magnitude  of  acceleration  that  can  be  prop¬ 
agated  beyond  it.  At  the  same  time,  the  critical  acceleration  always  de¬ 
creases  with  depth,  in  a  homogeneous  section  with  constant  slopes. 

27.  A  note  of  caution  is  in  order  for  dams  with  abrupt  changes  in  sec¬ 
tion  or  zoning  that  would  cause  a  reduction  in  yield  accelerations  for  slip 
surfaces  above  the  base  of  the  embankment.  For  example,  some  dams  ha-e  slopes 
that  steepen  abruptly  near  the  crest.  However,  for  upstream  slip  surfaces, 
the  reduction  in  yield  acceleration  due  to  steeper  slopes  is  usually  more 
than  offset  by  an  increase  due  to  lower  pore  pressures  if  the  steeper 
section  lies  above  the  pool  elevation.  Upstream  or  downstream  berms  will  also 
result  in  relatively  reduced  yield  accelerations  for  slip  surfaces  that  lie 
entirely  above  the  berms.  For  these  or  similar  cases,  a  profile  of  yield  ac¬ 
celerations  can  be  developed  from  stability  analyses;  Figure  11  can  be  used  to 
estimate  amplification  factors;  and  potential  displacement  can  be  calculated 
from  Figure  7  for  each  of  the  potential  slip  surfaces  identified  in  the  sta¬ 
bility  analyses. 

28.  The  authors  conclude  that,  except  for  a  few  special  cases,  deep 
sliding  surfaces  are  of  greatest  significance  when  evaluating  the  possibility 
of  displacements  that  could  threaten  the  integrity  of  the  structure,  and 
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accordingly,  they  look  for  a  limiting  amplification  factor  representing  slid¬ 
ing  surfaces  at  the  base  of  the  embankment.  Figure  11  shows  that  the  value  at 
the  mean  plus  2a  limit  is  approximately  3.0.  Applying  this  amplification  fac¬ 
tor  to  the  N/A  value,  0.17,  which  gives  an  upper  bound  of  1-m  displacement 
for  the  rigid-plastic  sliding  block,  the  ratio  of  critical  acceleration  to 
peak  bedrock  acceleration  is  0.5. 


PART  III:  CONCLUSIONS 


29.  The  results  of  analysis  of  earthquake  strong-motion  records  using  a 
sliding  block  model  and  a  decoupled  elastic  response  analysis  show  that  per¬ 
manent  displacements  for  deep-seated  sliding  surfaces  limited  to  less  than  1  m 
can  be  assured  if  the  ratio  of  critical  acceleration  to  peak  bedrock  accelera¬ 
tion  is  at  least  0.5.  This  value  is  considered  to  be  very  conservative  and 
subject  to  downward  revision  as  better  understanding  of  elastic-plastic  ampli¬ 
fication  response  of  embankments  is  developed. 

30.  Furthermore,  a  pseudostatic,  seismic  coefficient  analysis  can  serve 
as  a  useful  screening  procedure  to  separate  dams  that  are  clearly  safe  against 
earthquake-induced  sliding  failure  from  those  that  require  further  analysis. 
The  permanent  displacement  analyses  described  in  this  report  provide  a  ra¬ 
tional  basis  for  choosing  the  value  of  the  seismic  coefficient. 

31.  The  suggested  procedure  is  as  follows: 

a.  Carry  out  a  conventional  pseudostitic  stability  analysis,  using 
a  seismic  coefficient  equal  to  one-half  the  predicted  peak  bed¬ 
rock  acceleration. 

b.  Use  a  composite  S-R  strength  envelope  for  pervious  soils,  and 
the  R  (undrained)  strength  for  clays,  multiplying  the  shear 
strength  in  either  case  by  0.8. 

c.  Use  a  minimum  factor  of  safety  of  1.0. 

32.  This  procedure  should  not  be  used  in  the  following  cases: 

a.  Where  areas  are  subject  to  great  earthquakes  (of  magnitude  8.0 
or  greater) . 

b.  Where  materials  in  either  the  embankment  or  foundation  are  sus¬ 
ceptible  to  liquefaction  under  the  design  cyclic  loading. 

c .  Where  the  available  freeboard  is  small,  or  where  the  dam  has 
safety-related  features  that  art  vulnerable  to  small 
deformations . 
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Jennings,  Housner,  and  Tsai  (1968). 

Seed  and  Idrlss  (1969). 

Values  in  parentheses  are  for  reversed  direction  of  shaking. 


